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Abstract 

We study the propagation of helicity-1 gravitons in the Fierz-Pauh massive gravity 
in nearly Minkowski backgrounds. We show that, genericahy, there exist backgrounds 
consistent with field equations, in which the propagation is super luminal. The relevant 
distances are much longer than the ultraviolet cutoff length inherent in the Fierz-Pauli 
gravity, so superluminality occurs within the domain of validity of the effective low 
energy theory. There remains a possibility that one may get rid of this property by 
imposing fine tuning relations between the coefficients in the non-linear generalization 
of the Fierz-Pauli mass term, order by order in non-linearity; however, these relations 
are not protected by any obvious symmetry. Thus, among others, superluminality is a 
problematic property to worry about when attempting to construct infrared modifica- 
tions of General Relativity. 

1 Introduction and summary 

Lorentz-invariant massive gravity in four dimensions — the Fierz-Pauli theory — is of 
interest from the viewpoint of understanding the problems that may arise when one at- 
tempts to modify General Relativity in the infrared domain. Indeed, at the linearized level 
about Minkowski background, classical Fierz-Pauli theory exhibits the van Dam-Veltman- 
Zakharov discontinuity [I112] and non-linearity at large distances from gravitating bodies [3], 
whereas at the quantum level this theory becomes strongly coupled at unacceptably low 
energies [1]. In slightly curved backgrounds, the Fierz-Pauli gravity and its non- linear gen- 
eralizations necessarily contain an extra, Boulware-Deser mode [5], over and beyond the five 
modes of massive graviton; furhermore, one of the modes is always a ghost. These properties 
of the Fierz-Pauli massive gravity are reviewed, e.g., in Ref. [6]. 
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In this note we point out that the Fierz-Pauh theory has yet another pecuhar feature. 
Namely, we find that, generically, there exist shghly curved backgrounds consistent with 
field equations, in which some of helicity-1 graviton modes are superluminal. We recall that 
these modes are not pathological otherwise: they do not exhibit the vDVZ discontinuity, 
Vainshtein non-linearity or particularly bad UV behavior at the quantum leve£]. We check 
that the superluminal propagation occurs over large distances as compared to the UV scale 
intrinsic in the Fierz-Pauli gravity, so this phenomenon shows up within the domain of 
applicability of the effective low energy theory. 

In this regard, the Fierz-Pauli massive gravity is analogous to the Dvali-Gabadadze- 
Porrati model [7], which also has a superluminal mode in some legitimate backgrounds [5]. 
As discussed in Ref. [8], superluminal propagation signals that there is no UV completion 
of an effective low energy theory into quantum field theory or perturbative string theory, so 
the property we have observed is yet another feature to worry about. 

Our analysis has a loophole, though. There remains a possibility that the superluminal 
propagation of helicity-1 modes may be avoided by imposing an infinite set of fine tuning 
relations on the parameters in the non-linear generalization of the Fierz-Pauli mass term, 
order by order in the degree of non-linearity. These relations, however, are not protected by 
any obvious symmetry, so we do not find this "solution" attractive. 

It is worth stressing that we study here massive gravities that have the Minkowski metric 
as a solution to the field equations. The situation may be less pathological in theories whose 
"natural" background is different from Minkowski; a well studied example is massive gravity 
about (anti) de Sitter space-time [9l [TOl [IH [HI [131 [H] . 

Once Minkowski space-time is a solution to the field equations, the general mass term is 
a polynomial in {g^i, — r]f^u)- The original Fierz-Pauli action is quadratic, 

Sfp = ^ I d'x {-r^^^r^'^P^g^, - ^^,){g,^ - ^,^) + [^^'^{g^, - r^^,)]2} . (1) 

Our purpose is to evaluate the speed of helicity-1 modes of metric perturbations in nearly 
Minkowski backgrounds. So, we will make use of the perturbation theory in (^^^, — r]fj,^), 
where 'g^^^ is the background metric. To this end, we will have to generalize the above 
expression by adding cubic, and then higher-order terms. It is at this stage that the fine 
tuning "solution" of the superluminality problem shows up. 

We begin with studying in section [2] the Fierz-Pauli theory and its generalizations in 
cosmological backgrounds of small space-time curvature. In these backgrounds, the metric 
perturbations decompose into the tensor, vector and scalar modes with respect to spatial 
rotations. To set the stage, we consider in section 12.11 a particular form of the mass term 

"'^Helicity-0 (longitudinal) modes may be superluminal too, but they are less interesting in our context 
because of the other pathologies inherent in the longitudinal sector. 
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([T]) and find that the vector perturbations in this case are superluminal in fairly general 
backgrounds. By itself, this is not a particularly strong result, however. The degree of 
superluminality (c^ — 1), where c is the speed of the propagation of the vector modes of 
metric perturbations, is proportional to the deviation of the background from Minkowski 
space-time, {g^u — Tj^^). This implies that the contributions to (c^ — 1) due to O [{g^u — Vfj.u)'^]- 
terms in the action are of the same order. We consider the cubic terms in section 12.21 and 
find that to the first order in (^^^, — rjf^,^), the total (c^ — 1) may be made equal to zero by fine 
tuning a parameter in the cubic action. To see that similar fine tuning is required in higher 
orders in [g^^ — rj^^), we have to study the action containing higher orders in {g^^ — rj^i,). 

Proceeding with cosmological backgrounds is technically challenging at this stage. In- 
stead, we make use of the Stiickelberg formalism. We consider backgrounds of a special form 
and gravitons traveling along a particular direction. This enables us to separate helicity-1 
modes from other metric perturbations in a consistent way. In section 13.11 we cross check 
by rederiving, within this formalism, the fine tuning relation ensuring that = 1 at the 
linear order in {g^^^ — r]^^). In section [3^ we show explicitly that for the backgrounds we 
consider, the requirement that helicity-1 modes do not propagate in superluminal way gives 
rise, at the quadratic order in {g^^ — rj^y), to further fine tuning relations, now involving the 
coefficients in the fourth-order action. It is then straightforward to convince oneself that fine 
tuning proliferates to higher orders. We have found no arguments suggesting that our fine 
tuning relations are sufficient to avoid superluminal propagation of helicity-1 modes in arbi- 
trary backgrounds; neither have we found backgrounds that would rule out the fine-tuning 
"solution" to the problem of superluminality in the helicity-1 sector. 

2 Cosmological backgrounds 

The theory we consider in this paper is the Fierz-Pauli model whose action is 

S = SqR + Sm 

where Sm is the mass term including the quadratic part ([1]) and possible higher order terms, 
while Sgr is the action of General Relativity. In the further analysis we set IQttG = 1. 
Our purpose in this section is to study the propagation of gravitons in the spatially fiat 
cosmological backgrounds, in the regime when spatial momenta and frequencies are much 
higher than the graviton mass. In the coordinate frame where rj^j^^ in the mass term is 
diag[l, —1, —1, —1], the spatially fiat FRW metric is 

ds^ = n^{t)a^{t)dt^ - a\t)dx^ (2) 

In what follows, it is convenient to use conformal time t] related to time t by 

drj = n(t)dt 
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In the coordinates (77, x), the speed of hght is equal to 1. We will consider nearly flat metric, 
for which Sn = n — 1 and Sa = a — 1 are small. 

The field equations for the background, written in terms of conformal time, read 

f-V=H^ = eo, 2--n' = es (3) 
\a J a 

where prime denotes d/dr], the contributions eo , are due to the mass term and are of 
order 

eo , e, = ml [0{a - 1) + 0{n - 1)] 
The lowest order terms come from the quadratic part (JT]); their explicit expressions are 

6o = -^(a^-l), = -^(2(a^ - 1) + (aV - 1)) (4) 
2 2n 

where the expansion up to the first order in 5a and 5n is understood. 

Now, we can freely choose 5n and 5a at a given moment of time. Then these variables 
will change in time at low rate determined by the value of mc, which we take small. We will 
consider the time scales shorter than m'^^ so we treat 5n and 5a as constants. It is worth 
noting that the consistency of the system ([3]), i.e., covariant conservation of the effective 
energy-momentum tensor coming from the mass terms, implies that n{ri) obeys the equation 

-^n = -—a +n{es-eo) . (5) 

However, this equation determines the evolution of 5n and does not prohibit one to choose 
arbitrarily the value of 5n at a given moment of time. So, for our purposes we can indeed 
treat 5n and 5a as independent free parameters of the background. 

Let us now include metric perturbations, and write the perturbed metric in the following 
form, 

ds"^ = a^{t) [n^{t){l + hoQ)df + 2hQin{t)dtdx' + {-5ik + hik)dx'dx''] . 

Due to the 0(3) symmetry of the FRW metric one can study different helicity sectors sep- 
arately. We will be interested in vector (helicity-1) perturbations, for which, using the 
standard notation, we write 

hij = dWj + djWi, diW, = 0, 
hm = 0, hoi = Ui, diUi = 

The question is whether or not the propagating vector modes are superluminal for some 
choice of 5a and 5n. 
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2.1 Quadratic mass term only 

To warm up, let us consider the particular form of the graviton mass term given by ([1]), with- 
out cubic term. The complete quadratic action for vector perturbations in the background 
obeying the equations of motion ([3D is 

2U,AW- - UAUi + W-m^i + 2nW[m^i + Qe^UiU, + 2esWiAWi 

where A = didi. The last term in the above expression comes directly from the mass term 
in the action. In the case ([T]) the explicit expressions for Sm^, 6M'^ to the first order in 6n, 
6a are 

6m^ = ml{l + A6a + 26n), 6M^ = ml{l + A5a) 
It is convenient to define 

= 4e, + -^6M\ = 12eo + -^6m^ 

Then the linear equations for metric perturbations are 

M^AW - AHAU - 2AU' + 2AW" + AHAW = 

m'^U + 2AW' - 2AU = 

We now concentrate on frequency and momenta exceeding Hq and rriG, and make the Fourier 
transformation ^ iuo, A — p^. In this limit we obtain the following dispersion relation 
for the vector perturbations, 

To the first order in 6n, 6a we find 

This implies that for any background with 6n < 0, the speed of vector perturbations exceeds 
the speed of light. 

Importantly, — 1 as given by ([7j) is of the zeroth order in rriQ. On the other hand, 6n' 
is of the first order in this parameter. So, at small there is enough time for the vector 
perturbations to propagate, in superluminal way, before the space-time metric changes. In 
detail, the metric does not change during the time interval r ~ m^^ ■ |5a|^/^, where the factor 
|(5a|^/^ is obtained by inspecting eqs. (jl]) and ([5]). The distance the vector perturbations 
advance light in this time interval is L = (5c ■ r ~ m^^ ■ |(5a|^/^|(5?2|, which for not very small 



(2) 
tot 



E.n.M. 



a^drfd^x 



2 



m 



- 1 = -46n 



(6) 



(7) 
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6a and 5n is much greater than the strong couphng scale [1] inherent in the Fierz-Pauh 
theory, A^y = (m^Mp;)~^/^. Hence, the superluminal propagation is not hidden below 
the UV distance scale. This reasoning applies word by word to massive gravities with more 
general actions, to be discussed below. The effect of the superluminal propagation we analyze 
in this paper does occur within the domain of applicability of the effective field theory. 

2.2 The cubic action 

To the linear order in 5n, 5a, the speed of vector perturbations may receive contribution from 
cubic in {g^j.^ — rj^i,) part of the mass term. One may thus wonder whether superluminality 
may be eliminated by the proper choice of the parameters in this part of the action. To 
examine this issue, let us calculate the linear in 6a, 6n contribution to for general cubic 
term. As before, the mass term in the action is assumed to be Lorentz-invariant, contain no 
derivatives and have fiat metric as a solution to the massive gravity equations. Its quadratic 
part is of the Fierz-Pauli structure. The general expression for the cubic part is 

= J d\ {C-Al + 8- A^A^ + ^ ■ A3) , (8) 

where 

with C, S, T being arbitrary coefficients. The total mass part of the action is the Fierz-Pauli 
term ([T]) plus S'^ . 

The expressions (jS]) - ([6]) are still valid, but bm and bM are to be recalculated. The 
linearized expressions are 

W = m^(l + 45a + 25r2)-8£(45a + H-12J^(25a + H 
= m^(l + 4(5a) -8^(4(5a + 5n) -24J^5a 

The previous result ([7]) then modifies to 

c^-l= ^-^~^g .45n (10) 

Thus, to avoid superluminality at the first order in (5a, 8n, one has to impose the fine tuning 
relation 

3jr = ^2^ (11) 

We will reproduce this fine tuning relation in section 13.11 in the Stiickelberg approach, and 
in section 13.21 we proceed to higher order terms. 
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3 The Stiickelberg approach 



3.1 Cubic order fine tuning revisited 

To continue the analysis, it is convenient to use the Stiickelberg formalism, wich will also 
allow us to study backgrounds other than cosmological. Let us write the metric as follows, 

gfiiy = Vi^u + hf,u + hf^u , (12) 

where the first term is the Minkowski metric, the second term corresponds to non-trivial 
background and is assumed to be small, and the third one describes perturbations about 
this background. 

To perform the Stiickelberg analysis, one enlarges the set of fields in the theory by 
introducing new fields S,'^ and g^^ in the way dictated by the gauge symmetry of General 
Relativity [H [151 [T6|, [T7], see also Ref. [6] for detailed discussion, 

gA^) = gf.u{x + + ^M^^ 9ux{x + + g,^xix + + d^^^ hpi^ + . 

We are going to consider high graviton momenta and slowly varying backgrounds, so we 
neglect the derivatives of h^,^. Once the field g^u{,x) is gauge fixed, mixing between its 
fluctuations and the field becomes irrelevant at high enough momenta, and we are left 
with the theory of one vector field determining the interesting part of metric perturbations 
via 

Ku = d^i, + d,i, + d^ehx, + d^ehx, + d^ed.ix + Kd.edui' ■ (13) 

Throughout this paper indices are raised and lowered with the Minkowski metric. Note that, 
nevertheless, this expression is exact for slowly varying backgrounds, i.e., no expansion in 
h^y has been made yet. 

The field does not enter the Einstein-Hilbert part of the action, and obtains kinetic 
term from the mass term in the original action. To evaluate its speed, we study the action 
to the quadratic order in S,^. We begin with the linear order in h^y\ in this way we are going 
to reproduce the fine tuning relation fill I) in the Stiickelberg formalism. Hence, we need the 
mass term to the cubic order in {g^^ — rj^y), 

^^2+3) _ ^ g{3) ^ j ^4^ . (_^2 ^ ^ Q ■ Af + S ■ A^A^ + ^ • A3] , 

where 

A ^ -mil A , (14) 

and the combinations Ai, A2, A3 are given by (111). According to the above discussion, we 
will need the terms in this action that are quadratic in the Stiickelberg field and zero and 
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first order in h^,^. Note that these terms appear both due to the quadratic and cubic terms 
in the action and due to the non-hnearity of the gauge transformation (fT3|) . 
Making use of and ([13]) one obtains 5^+^^ = / d'^x with 

= {-2^ [{d.i^f + - {2A + 3^) ijr^'d.d^e 

+2{A£ + 3^ - 2A) dx^^ ■ h^'d^i, 

+2{A - E) K,Ui^ + 2{£ + 6C) h{d^ef - (3^ + 4^) V^^nr} (15) 

where 

h = h^^r]'''' . 

Let us consider the frame in which the Minkowski metric entering the mass terms has the 
standard form rj^u = diag[l, —1, —1, —1] and study the backgrounds whose space-time metric 
is diagonal in this frame, h^^i, = diag[/ioo; hu, /122, ^33], so that 

ds^ = {l + hoo)dt^ + + h^i){dx'y . (16) 

i 

For time-independent hf^u this metric is a solution to the Einstein equations, so in massive 
theory the derivatives of h^^^, are proportional to ttiq. Hence, it is legitimate to consider this 
metric as slowly varying. This justifies the use of the expression f[T5]) . and also enables us to 
neglect the derivatives of /i^j^ when solving the equations for 

Let us further specify to gravitons propagating along the 1-st axis, whose momenta are 



Po = Pi= pS 



il 



In this case the action (fT5|) is a sum of three independent parts: the action that involves the 
field only, the action for and the action for the pair (^°, ^^). An obvious reason for this 
decoupling is that the matrix h^^i, does not mix, say, with other components of while 
d^C,'^ involves and only. Hence, the transverse polarizations ^'^ and decouple from 
each other and from longitudinal polarizations. These are precisely the (linear combinations 
of) helicity ±1 modes we are interested in. Their dispersion relations are readily calculated. 
One finds that the dispersion relation for both and to the linear order in /i^^, is 



— — — (/loo + rill J 



The physical speed of these modes (recall that they propagate along the 1-st axis) is 



where 

• _ 900 _ 1 + hpQ 
\9ii\ l-hu 

Hence to the hnear order in h^^^, the physical speed of hehcity-1 gravitons is given by 

c = 1 - ^ — ^ — -{hoo + hi) 

where {hoo + ^ii) niay have either sign. Recalhng (fT4|) and having in mind the two forms 
of background metric ([2]) and (fT6|l . we see that the expression (fT9l) coincides with (fTOj) . In 
this way we recover the resuh of section 12.21 the condition for the absence of superluminal 
propagation of hehcity-1 modes at the hnear order in /i^j, is 3JF = —AA = tOq. 

3.2 Fine tuning at fourth order and beyond 

Let us now study the theory tuned according to ( ITTi) . At the first order in h^u, hehcity-1 
gravitons in this theory propagate at high enough momenta precisely with the speed of light. 
We will now see that at the second order, superluminlity of helicity-1 gravitons generically 
reappears, unless one imposes additional fine tuning relations. 

Let us introduce the general fourth-order mass terms and expand the action up to the 
order oi 0{h x h"^). The fourth order terms are 

S^^^ = j d'x [X-{A,) + J-iAiY + }C-iAiAs) + C-iAlA2)+M{A2f] . 

Here, in analogy to ^2,3 in ([8]), the irreducible term is 

A4 = ry^2!^i _ _ _ ■^'^i^'i . (^g^^^^ — r^^^yj . . . ((7^41^4 — rjfi^u^) 

and I, . . .Ad are arbitrary constants. Plugging the Stiickelberg decomposition f|T3l) into this 
action, and keeping the terms of order 0{h x h) in the complete action, we obtain the 
contribution to the quadratic action of the Stiickelberg fields in the form Sm^ = J d^x C^'^\ 
where 

(h, = - (3C + 2C) (h'^n^,) -iS + AM) {K^V' ■ ^0^^) - {2£ + 3/C) {Wh"'d,d,iy) 

- (3^ + 4T) {eh^'h%d,ix) - {2A + 6^ + 4X) {h^^V^e^iu) 

+ 2(3^ + 4f + 4J + 6/C) {d^i^ ■ VX'd^ip) + 2(12:r + C) h\d^ef 

- (3/C + 4f - 2A) (hh^^COC) + (6/C + 16£ + 24C + AS) (hdx^^ ■ h^'^d^^^) 

+ 2{-A + 3J^ + A£ + 8M + 21) (F'S^e.)' - i2A + QJ^ + 41) ieh''%xd^d,^^) 

The total Lagrangian is the sum of this term and (ITSl) . We continue to consider the diagonal 
background metric (fT6!) and gravitons propagating along the 1-st axis. The fields and C,^ 

9 



(18) 
(19) 



again decouple, for the same reason as in section [STTl Making use of (ITTl) . we obtain the 
following dispersion relation for the coordinate frequency of the mode to the second order 
in h,. 



UJ^ = {l + {hoo + hn) 

- 9T_ 9 /I - ZL>^ - 'Ar _ A - T _ _1_ _ 

(/loo + ^ll) f - P ■ 



+ 



- 3A-2I- 2A-AS-3JC- ^A-I ,^ ^ , 
A — 2A A ~ 



Finally, expanding n ^, defined according to (ITSl) . to the second order in h^i, we find that 
the physical speed ( JTTll of the helicity-1 graviton with the polarization is given by 



3A-2I- 2A-A£-3IC^ 3A-2X - ' 
— ~A — ^ 2A ^ — ~A — ^ ^ '* 



(/loo + /ill) 



Since the quantities (/loo + ^ii)-, (/'-oo ~ ^ii); ^22 and /i are independent of each other, we 
see that to avoid the superluminal propagation, one has to impose additional fine tuning 
relations 

3^ - 2T = , 2^ - 4£ - 3/C = . (20) 

The same relations ensure that the mode ^ is not superluminal in the background we 
consider. 

The analysis in this section reveals the following property: the dispersion relation for 
helicity-1 gravitons involves the coefficients of the irreducible higher order terms, which are 
multiplied by sign-indefinite combinations of the background metric. At the third order, 
the relevant coefficient is JF, the irreducible term is A3, and the combination of metric is 
(/loo + ^11); while at the fourth order, these are X, A4, and /i22(/ioo + ^ii); respectively. 
Unless these coefficients are fine tuned, helicity-1 gravitons propagate in superluminal way 
in backgrounds with appropriate signs of the background metric coefficients /i^j^. Indeed, 
the relation fllip contains JF, and the first of the relations (I2U]) contains X. One can check 
that this property holds at higher orders. Hence, the superluminal propagation reappears at 
higher orders unless at least one fine tuning relation is imposed at each order in ((7^,^ — ?7^j,) 
(in fact, our fourth-order analysis shows that the number of fine tuning relations is larger 
than one). 

There is no obvious symmetry behind the relations like ( fTTj) and ( l20l) . Thus, helicity-1 
gravitons propagate superluminally in the Fierz-Pauli massive gravity unless this theory is 
heavily fine tuned. 

This work has been supported in part by Russian Foundation for Basic Research, grant 
08-02-00473. 
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